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Abstract. We report on full-polarization radio observations of the Perseus cluster (Abell 426) using the Westerbork Synthesis 
Radio Telescope (WSRT) at wavelengths from 81-95 cm. We detect faint, very extended polarized emission throughout the 
cluster region. We have employed a novel technique. Rotation Measure synthesis iBrentiens & de Bruvn 2005) to unravel the 
polarization properties of the emission across the full field of view. We detect polarized emission over a wide range of RM from 
about to 90 rad nT^. Low RM emission (RM < 15 rad m"-) is attributed to the local Galactic foreground. It has a chaotic 
structur e with smooth changes in polarization angle on scales of the order of 10' -30', not unlike those seen bv lHaverkorn et all 
j2003ah at the same frequencies. Emission at values of RM > 30 rad m ^ on the other hand, shows organized structures on scales 
up to a degree and displays rapidly fluctuating polarization angles on scales of the synthesized beam. A Galactic foreground 
interpretation for the high RM emission can not be ruled out, but appears extremely implausible. WSRT observations at 21 cm 
of the RM of a dozen discrete sources surrounding the Perseus cluster indicate a smooth large-scale gradient in the Galactic 
foreground RM. The diffuse structures have a clear excess RM of about 40 rad m"- relative to these distant radio galaxies. 
This excess Faraday depth, the generally good spatial association with the cluster and the different morphology of the high RM 
emission, compared to the genuine Galactic foreground emission, all point to an association of the high RM emission with the 
Perseus cluster. The polarized emission reaches typical surface brightness levels of 0.5-1 mjy per 2' x 3' beam and must be 
rather highly polarized 20%). Due to dynamic range limitations and lack of sensitivity to large-scale structure we have not 
yet detected the corresponding total intensity. Most of the polarized emission, located at distances of about 1° from the cluster 
centre, appears too bright, by about 1-2 orders of magnitude, to be explainable as Thomson scattered emission of the central 
radio source off the thermal electrons in the cluster. However, this remains a viable explanation for the highly polarized 2 1 cm 
emission from the inner 10'-20'and part of the 81-95 cm emission. The bulk of the emission associated with the Perseus cluster 
may instead be related to buoyant bubbles of relativistic plasma, probably relics from still active or now dormant AGN within 
the cluster. A lenticular shaped s tructure, referred to as the lens, and measuring 0.5-1 Mpc is strikingly similar to the structures 
predicted by EnBlin et al. 1 1998). At the western edge of the cluster, we detect very long, linear structures that may be related 
to shocks caused by infall of gas into the Perseus cluster along the Perseus-Pisces filamentary structure of the cosmic web. 

Key words. Galaxies: clusters: general - Galaxies: active - Polarization - Magnetic fields - Techniques: image processing - 
Radio continuum: general 



1. Introduction 

Clusters of galaxies play a crucial role in many fields of extra- 
galactic research. Located at the crossroads of the filamentary 
structures known as the cosmic web, they point to the areas 
with the deepest potential wells as predicted in the standard 
CDM model for structure formation in the Universe. In addition 
to the dominant dark matter and the large numbers of visible 
galaxies, clusters also contain enormous reservoirs of diffuse 
hot gas. This gas is visible directly through their X-ray emis- 
sion, and indirectly through their dynamical effects on galaxies, 
such as stripping, and the shaping of head-tail radio sources. 

Because of their deep potential wells, clusters continue to 
accrete gas from their surroundings. At the cluster-IGM in- 
terface this should lead to (mild) shocks CBurns.,1998J which 
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might be observed through enhanced X-ray and radio emission. 
In the cluster Abell 3667 infall indeed appears to be responsi- 
ble for the generatio n of very large radio and polarized source 
structures (Rottgerin g et al.ll9 97: Johnston-Hollitt 2004). 

The large-scale distribution of the hot gas can in principle 
be derived from the X-ray surface brightness profiles. However, 
there are two caveats. First, due to the quadratic dependence of 
X-ray emission on gas density, this distribution depends on the 
level of dumpiness in the medium. Furthermore, temperatures 
below a few times 10^ K also lead to reduced X-ray emission 
and therefore might lead to an underestimate of the amount 
of gas. Recent XMM observations ( Kaastra et al. 2003) indeed 
suggest that significant amounts of lower temperature gas ex- 
ist at the outskirts of clusters. These observations may have a 
direct bearing on the question 'Where are the missing baryons 
in the universe?' (iCen & Ostriker..l999.) . An independent Un- 
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ear probe of the gas density would be helpful to unambigu- 
ously deduce the gas distribution. The scattering of the 2.7 K 
microwave background photons by the cluster electrons, the 
Syunyaev-Zeldovich effect, is such a method. Another linear 
probe of the gas density is provided through Tho mson scat- 
tering of radio emission loc ated within the cluster dSvunvaevI 
Il982t IWise & Sarazinlll990l) . This method, which has a fair 
number of assumptions built in, has a better chance of detec- 
tion in low redshifts clusters because of confusion problems. 
Until now no unambiguous identification of Thomson scatter- 
ing in clusters has been made. Very high d ynamic range 2 1 cm 
continuum observations with the WSRT (de Bruyn 1995) may 
have detected such emission and were the original motivation 
for the observations described in this paper. 

Recent radio studies of clusters have drawn attention to 
another interesting aspect of clusters, namely the role played 
by radio loud AGN in the thermal balance of the intracluster 
medium. Combined radio and X-ray images of the inner parts 
of the Perseus and other rich clusters w ith strong radio sources 
llFabiaii et alJl20Q3[ Iciarke et alJl20Q4l) clearly show how the 
high pressure radio lobes of the central radio sources are the 
cause of cavities in the X-ray surface brightness. 

If AGN activi t y in the radio band is episodic 
dSchoenmakers et all 2000l) . clusters should contain many 
relic, or fossil, radio source s. Such sources have now been 
detected in several clusters llF.nRlin et al.lll998t iGovoni et al.l 
l200l'.'20Q5').'EnBl in & Gopal-Krishnal ( l200ll) show that orphan 
radio relics would slowly rise in the cluster atmosphere as a 
result of buoyancy forces. The adiabatic energy losses (on both 
particles and magnetic field) would then render these relics 
efl'ectively invisible at GHz frequencies. These 'bubbles' could 
then 'hibernate' for periods of several Gyrs and be invisible 
except at extremely low frequencies. These relics could be 
'woken up' by the shocks associated with the inflow of gas 
from the surroundings. 

In this paper we describe observations of polarized radio 
emission from the Perseus cluster which provide a new tool to 
study the three components of the diffuse intracluster medium: 
thermal gas, relativistic gas, and magnetic fields. Low fre- 
quency radio polarization data may convey information about 
extremely low density regions if they are pervaded by a mag- 
netic field. Such measurements, however, are complicated by 
instrumental and astrophysical effects which have limited their 
application in the past. First, off-axis instrumental polarization 
of radio telescopes is generally quite severe at low frequen- 
cies. Using mosaicking techniques to image large fields will 
alleviate this to a large degree. Second, ionospheric Faraday 
rotation can rotate polarization angles by up to hundreds of de- 
grees during the synthesis time at a frequency around 325 MHz. 
This can be limited by observing at night and using GPS-based 
and ionosonde-based estimates of the ionospheric electron con- 
tent. Third, the percentage polarization of radio sources is gen- 
erally low at low frequencies as a result of internal depolar- 
ization due to Faraday rotation within the emission regions 
jStrom & Conwav 1985: Conway & Strom 1985) and as a re- 
sult of beam depolarization. When obse rving sourc es with a 
large rotation measure, or Faraday depth ('Rurn ll96fil) . a fourth 
comphcation arises: bandwidth depolarization. 



In order for standard polarization analysis (e.g. RM-fitting) 
to be meaningful a minimum S/N of about 3 is required in 
the linearly polarized signals for each frequency. Increasing the 
bandwidth to improve the S/N, however, conflicts with the re- 
quirement to limit the bandwidth to avoid bandwidth depolar- 
ization. For example, observing a polarized source at 350 MHz 
with a bandwidth of 80 MHz, limits the allowed RM to about 
5 rad m - before Faraday rotation quenches all polarized sig- 
nals by an order of magnitude. To obtain good sensitivity yet 
avoid bandwidth depolarization requires a wideband correla- 
tor with a large number of channels. As of April 2002 the 
Westerbork Synthesis Radio Telescope (WSRT) can cross cor- 
relate, in full polarization, data from all 14 telescopes in up 
to 1024 spectral channels over an 80 MHz band. Employing 
a tec hnique, which we have cal led Rotation Measure synthe- 
sis ( B rentiens & de Bruvnl2005l) and which is briefly described 
below, then provides an elegant and powerful method to allow 
noise limited polarimetry using a wide frequency band with 
arbitrarily low S/N ratio in each individual narrow frequency 
channel. 

A powerful feature of RM-synthesis is the ability to search 
for weakly polarized emission over a wide fie ld of view. Asim- 
ple version of the method was developed by Ide Bruvr] ( 1 19961) 
for the old WSRT 8 channel continuum backend. It was imple- 
mented in the WSRT reduction package NEWSTAR and first 
applied to data fro m the highly polariz ed millisecond pulsar 
PSR J0218-H4232 (Nava rro etalJ ll995). The method will be 
described briefly in Section |31 but a full derivation, discussion 
and analysis of its properties is presented in a companion paper 
dBrentiens & de Bru vn''200.5Y 

The outline of this paper is as follows. In Section |2l we 
present the deep WSRT observations of the Perseus cluster mo- 
tivating this paper Section |3 describes the total intensity and 
polarization calibration. Sectionl^describes how the basic data 
products were made. Section|5ldiscusses instrumental artefacts 
and how they can be distinguished from the interesting astro- 
nomical signals. Sections |6l and discuss the nature and lo- 
cation of the various types of polarized structures. Section |S1 
presents our interpretation of the features that we attribute to 
the Perseus cluster. Some puzzles and planned future observa- 
tions are presented in Section|9lwhile SectionFTOIconcludes. 

2. Observations 

The observations were conducted with the WSRT. The array 
consists of fourteen 25 m dishes on an east-west baseline and 
uses earth rotation to fully synthesize the uv-plane. For full 
imaging over the whole primary beam it takes six array con- 
figurations in which the four movable telescopes are stepped at 
12 m increments (i.e. half the dish diameter) with the shortest 
spacing running from 36 m to 96 m. This provides continuous 
uv-coverage with interferometer baselines ranging from 36 to 
2760 m. The 12 m increment creates an elliptic grating lobe 
with a radius in Right Ascension of about 4° from the phase 
centre, which places it well beyond the 10 dB point of the pri- 
mary beam which measures 2°4 full width at half power. Self 
confusion is therefore not a problem within the -5 dB point of 
the primary beam. The angular resolution in Right Ascension is 
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Observation ID 


9-A (m) 


Stai-t date (UTC) 


End date (UTC) 


10208707 


72 


2002/12/02 16:06:20 


2002/12/03 04:04:50 


10208767 


36 


2002/12/06 15:50:30 


2002/12/07 03:49:00 


10208819 


96 


2002/12/10 15:34:50 


2002/12/11 03:33:20 


10208868 


84 


2002/12/12 15:27:00 


2002/12/13 03:25:30 


10208937 


48 


2002/12/14 15:19:10 


2002/12/15 03:17:40 


10208997 


60 


2002/12/17 15:07:20 


2002/12/18 03:05:50 



about 0.8'. The half power beam width, grating lobe and reso- 
lution all correspond to a frequency of 350 MHz, in the middle 
of the observing band, and scale inversely with frequency. The 
grating lobe and the angular resolution furthermore scale with 
l/sin(5) in the declination direction. Basic observational data 
are collected in Tabled The pointing and phase centre of the 
telescope was directed towards (B 1950.0): RA = 03''15™, Dec 
= +41°15', which is located between the three dominant radio 
sources in the Perseus cluster This position was chosen to be 
identical to tha t of previous WSRT observations described in 
ISiibrind ( ll993h to allow a comparison between the datasets. 

The feed/receiver is part of the new Multi Frequency 
FrontEnd package (http://www.astron.n1/wsrt/i and covers the 
frequency range from about 310-390 MHz. These frequencies 
correspond to wavelengths from 78-97 cm but for historical 
reasons, and brevity, we will continue to refer to this as the 
92 cm band. The 92 cm band is usually largely free of radio fre- 
quency interference (RFI) during evening, nighttime and week- 
end observing at the location of the WSRT. The full band can 
be completely covered by the new wide band correlator which 
can process 8 independently tunable bands of 10 MHz. Each 
band is covered by 64 channels in 4 cross-correlations to re- 
cover all Stokes parameters. We used Hamming tapering in the 
lag-to-frequency transform yielding an effective spectral reso- 
lution of 0.31 MHz. Hamming taper, rather than Hanning, wa s 
used to lower the distant spectral side lobe level llHarriJl978t) . 
The channel separation was 0.156 MHz. The frequencies of 
the 8 bands were centred at 319, 328, 337, 346, 355, 365, 374, 
and 383 MHz. A gap at 360 MHz was introduced to avoid lo- 
cal RFI. Gibbs ringing at the video edge of the band and the 
increased (digital) noise at the upper filtered end of the band 
led to useful data in 56 (out of 64) channels of each 10 MHz 
band (channels 3-58). Unfortunately, serious external RFI was 
encountered during several hours into the evening at the start 
of the observations. In several 10 MHz bands these intense but 
narrow RFI signals 'ring' through the whole band, despite our 
use of a Hamming taper. By combining pairs of odd-even chan- 
nels in the subsequent processing this ringing could be elimi- 
nated and the final analysis was therefore done using 28 chan- 
nel pairs for each 10 MHz band. Henceforth we will refer to 
such a channel pair as a channel. It has an effective spectral res- 
olution of about 0.4 MHz. The RFI frequently spoiled the total 
power monitoring data preventing a precise absolute calibration 
of the data. This absolute calibration was therefore done using 
the 1984 WSRT observations at 327 MHz (.Siibring. 1993,1 



Full syntheses with the WSRT usually need not be inter- 
rupted during a 12 hour track. Because of its equatorial mount 
- no change in parallactic angle - it also suffices to observe 
only one polarized and one unpolarized calibrator to calibrate 
instrumental leakages. Every observation was bracketed by two 
pairs of calibrators, one polarized and one unpolarized: 3C 345 
and 3C 48 were observed before the target, while 3C 147 and 
the bright, highly polarized eastern hot spot in DA 240 were 
observed at the end. The calibrator observations lasted 30 min. 
The time resolution in the data was 30 s, which was generally 
sufficient to fully sample ionospheric phase fluctuations with- 
out serious phase decorrelation effects. This time sampling was 
also sufficient to avoid time smearing for sources at the outer 
edge of the field. 



The system temperature of the WSRT at 350 MHz towards 
the Perseus cluster is about 125-150 K. With a net integration 
time of 72 hours this should result in a thermal noise of about 
20 yuJy for an effective bandwidth of 8 x 8.7 MHz. However, we 
did not reach this level for the following reasons. Calibration 
problems did not allow us to include the upper two frequency 
bands (374 and 383 MHz). Due to serious malfunctioning in 
two telescopes/receivers (RTO and RTC), one of which (RTC) 
provides half of the long baselines in the array, we decided to 
use only the inner half of the array in the final analysis. This 
also helped to render the very faint but very extended polarized 
emission visible. To avoid grating lobes, and the inherent self 
confusion once emission fills the whole primary beam, we did 
not use natural weighting in the imaging (i.e. the redundant 
baselines were not included in the imaging process). Hence the 
final image cubes were made with only 18 of the maximum 91 
baselines for each 12 h synthesis. This resulted in a synthesized 
beam of 2 .0' x 3 .0' . Henceforth, when quoting flux densities, we 
will always refer to this beam size unless otherwise specified. 



The theoretically achievable noise level should increase by 
a factor 2.5 to 50 yuJy (this corresponds to 50 mK brightness 
temperature). The final noise level, after flagging about 25% of 
the data that were affected by faint RFI or backend problems, 
was about 70 //Jy per beam. This, however, is still much lower 
than the classical confusion noise of about 1 .5 mJy for this fre- 
quency and angular resolution and presents by far the deepest 
low frequency image ever made. 
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3. Data reduction 

3.1. Total intensity calibration 

The Perseus clust er was first o bserved with the WSRT at 
327 MHz in 1984 ( ISiibrinJl993t) . The cluster contains a total 
of about 50 Jy of diffuse emission, concentrated in NGC 1275 
(3C 84) and NGC 1265 (3C 83.1), and an unusually rich con- 
centration of bright 4C sources at the edge of the primary 
beam (4C 41.08, 4C 42.09, 4C 43.09). The complexity of the 
brightness distribution necessitated several non-standard steps 
in the processing which was done using the WSRT-tailored 
NEWSTAR reduction package. Initially the total intensity self- 
calibration of the data did not converge satisfactorily. This we 
attributed to the absence of a sufficiently accurate initial ampli- 
tude calibration (no system temperature corrections were ap- 
plied) in the starting model. This is a well-known problem in 
WSRT observations of very complex fields, and is fundamen- 
tally due to the 1 -dimensional nature of the WSRT array. By us- 
ing the redundancy constraints in the array ( Wieringa 1992) we 
were able to generate an improved starting image which con- 
verged much better after self-calibration. Obviously bad data 
were edited at the start; more sophisticated flagging was done 
in an iterative fashion on the basis of the selfcal residuals. The 
total intensity image of the Perseus cluster is shown in Fig.^ 
This image was made from 1 1 channels in band 3 (average fre- 
quency 348 MHz) and is shown at fuU resolution (0.9' x 1.3' 
beam). The noise level varies across the image from about 0.3- 
1.0 mJy beam ' with the higher values in the inner 2°-3°of the 
imaged area surrounding 3C 84 and 3C 83. IB. Compared to 
the peak intensity of 18 Jy on 3C 84 this represents a dynamic 
range of about 20000:1. The observed noise is only a factor 
1 .5-2 above the classical confusion noise level. 

High dynamic range imaging at low frequency usually is 
limited by the phase stability of the ionosphere. In general the 
ionospheric phase fluctuations were rather modest during the 
6 nights of the observation, with peak excursions on a 1 .4 km 
baseline of about 30° on timescales of 3-5 minutes. In sev- 
eral nights, however, we experienced periods of a few hours of 
much faster, temporally unresolved, scintillation related phase 
excursions. These phase variations point to small scale struc- 
ture in the ionosphere and led to imaging problems, because 
they invalidate the standard self-calibration assumption of po- 
sition invariant errors. These manifest themselves in the im- 
ages as spiky patterns surrounding bright sources at the edge 
of the observed field (Note that the WSRT instantaneous re- 
sponse has a fan beam response which rotates clockwise from 
position angle 90° to 270 ° during the 12 h synthesis time). 
The most visible ionospheric problems are those surrounding 
4C 43.09, which is located about 1°8 NNE from the phase cen- 
tre. They are also visible around 4C 41.08, located about 1°6 
east of the cluster centre. We also detect weak amplitude er- 
rors around these sources which are due to pointing problems. 
The results described in this paper are concerned only with the 
diffuse polarization in the inner 3°of the cluster and are not af- 
fected by pointing. We have recently started to experiment with 
a new iterative 'source-by-source' self-cahbration procedure. 



called peeling in the context of LORA^R' , to remove very bright 
sources with their own ionospheric phase and telescope com- 
plex gain solutions. Because this paper is mainly concerned 
with the polarized emission we will not elaborate on this pro- 
cedure here any further. 

3.2. Polarization calibration 

The WSRT telescopes are equipped with a pair of orthogonal 
linear feeds. All four cross correlations between the incident 
signal are formed. The polarization c alibration followed the 
pro cedures described in the papers by iHamaker et all lll996l) 
and lSault et al ]lll996f). The overall (on-axis) instrumental leak- 
ages, typically 1-2%, were calibrated using the unpolarized 
calibrator 3C 147. The off-axis polarization behaviour is much 
more complicated and is described in more detail in Section 

The instrumental polarization corrections were transferred 
to the polarized calibrator source DA 240 which is then used to 
phase align the two orthogonal linear polarizations. The east- 
ern hot spot of DA 240, has about 600 mJy of linearly polar- 
ized signal and its RM is about +4 rad m"^. The final polar- 
ization corrections were then transferred to the target source. 
The 'peeling' method alluded to above was also used to correct 
for the position and frequency dependent instrumental polariza- 
tion. This was done on a channel-by-channel basis. Doing this 
for 6x28=168 channels made this a very time consuming pro- 
cess. For each channel we performed a full self-calibration as 
well as polarization calibration using the same starting model 
adjusted to take case of the different primary beam attenua- 
tions at the different frequencies. The self-calibration model 
included (instrumental) polarization for the brightest discrete 
sources. The 'peeling' was done for 4 different sources (3C 84, 
3C 83.1 A + 3C 83. IB, 4C 43.09 and 4C 41.08) in succession. 
The full details of this procedure, which are not very relevant 
for the purpose of this paper which deals mainly with the po- 
larization results, will be described elsewhere. 

At the observed frequencies ionospheric Faraday rotation 
can sometimes be a problem, especially during day-time ob- 
serving. We do not have a strongly polarized source within the 
Perseus field on which we can 'self calibrate' the ionospheric 
Faraday rotation. We therefore made use of data from two near- 
simultaneous 12 h syntheses of NGC 891, observed for a dif- 
ferent project during the 2-week observing campaign on the 
Perseus cluster. T his field contains the highly polarized pulsar 
PSR J0218H-4232 jNavarro et alJll995l) . The pulsar's polariza- 
tion angle revealed a Faraday rotation of about 60° at the start 
decreasing to about zero after 4 hours and remaining small for 
the remaining 8h of the observation. The Faraday rotation oc- 
curred in the few hours before and just after sunset. Because 
the declination and epoch of the NGC 891 and Perseus obser- 
vations were very similar we adjusted for the 1 h RA-difference 
and de-applied a polarization angle rotation of 45° decreasing 
to zero during the first 3 hours of all six 12 h Perseus cluster 
syntheses. 



' a name suggested by Jan Noordam 
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Fig. 1. Total intensity image of the Perseus cluster, observed with the WSRT in 2002, using data from band 3. Contour levels are 
at -2, 2, 4, 8, 16, ... , 16384 mJy beam"'. 



4. From Q-U image cubes to RM-cubes 

For the final polarization analysis we eliminated 42 channels 
from the available 168 because they showed faint traces of RFI 
or dynamic range problems related to specific telescopes. The 
remaining 126 images were carried through for further process- 
ing in RM-synthesis. The average noise level in the Stokes Q 
and U polarization images of a single channel was rather uni- 
form and varied from about 0.7-0.8 mJy. 

The wide range of frequencies in the data set impUes that 
the primary beam width changes significantly across the band. 



At the mid-band frequency of 345 MHz the half power beam 
width (HPBW) measures about 1°22. The WSRT primary beam 
is well approximated by a cos(cvr)* function where c is a con- 
stant equal to about 0.064, v is the frequency in MHz and r the 
radius from the pointing centre in degrees. For the polarized 
structures discussed in this paper, which are detected out to a 
radius of about 1?5, the maximum reduction is about a factor 
3.5 at the highest frequency of 370 MHz used in the analysis. 
The differential attenuation between 315 and 370 MHz at this 
distance from the pointing centre is still only a modest factor 
1.4. We have not tried to correct the data for these differential 
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effects. In a future paper we will return to the spectral prop- 
erties of the polarized emission and, if we manage to recover 
them, the associated total intensity structures. 

4.1. Rotation-measure synthesis 

The basic idea behind rotation measure synthesis is that one 
derotates, for every pixel in each channel image, the Q-U vec- 
tors in order to compensate for a certain assumed rotation mea- 
sure. After derotation, the channel images are averaged. This 
procedure maximizes sensitivity to radiation at the assumed 
Faraday depth, because that emission is coherently added. All 
other emission will add only partly coherently, hence the sensi- 
tivity to emission not at the assumed Faraday depth is red uced. 
The es sence of this procedure was akeady mentioned bv lBurnI 
(Il966l). The procedure is not unknown in the pulsar community 
(Mit ra et al. 2003; Weisberg et al. 2004). Howeve r it is usually 
only applied to a si n gle source or Un e-of-sight jBower et alJ 
Il999t lKilleenlfl998i) . lde BruviJ ( Il996l) applied the method to 
every pixel in WSRT 350 MHz radio synthesis data of the field 
surrounding the galaxy NGC 891 and the highly polarized pul- 
sar PSR J0218+4232, but with only 8 frequencies available 
the RM side lobes were rather high. Following the comple- 
tion of the new broadband (8x10 or 20 MHz) 250000 chan- 
nel WSRT backend the method was expected to become much 
more powerful. A full d escription of RM-synthesis is presented 
in a companion paper (Bre ntiens & de BruvrJbOO^ . Here we 
summarize the essential features of the method as applied to 
the Perseus cluster data. 

The derotation of the multi channel complex polarization 
images can be performed very efficiently in terms of computer 
time. One first computes the complex polarization P - Q + ill 
and r nultiplies with a complex pha se factor to perform the rota- 
tion. iBrentiens & de BruvrJ ll2005 l) derive that a general inver- 
sion of the polarization as a function of wavelength squared is 
given by: 



XT W(/i2)P(^2)e-2i^('''-''o>di2 



where (p is the Faraday depth. 



(1) 



0.81 



Jthe 



eB dl. 



(2) 



and F{(p) is the emission as a function of Faraday depth, * 
denotes convolution, R{(p) the rotation-measure transfer func- 
tion (RMTF), Aq is the square of the wavelength to which 
aU vectors are derotated, and W{A^) the sensitivity as a func- 
tion of wavelength squared. This function is also known as 
the sampling function, or weight function. F{(p) is measured 
in Jy beam ' rmtr'. For extended sources that are discrete 
in 0, this corresponds to the polarized surface brightness. For 
sources that are extended with respect to 0, it is the polarized 
surface brightness per RMTF beam width. 

If the bandwidth of an individual channel is much less than 
the total bandwidth of the observation, we may approximate 
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Fig. 2. The RMTF of the observations. The FWHM is 15.2 
rad nT^. The real part is the response parallel to the polariza- 
tion vector and the imaginary part is the response orthogonal to 
it. Because of an uncertainty of the Faraday depth of the peak 
of typically a few rad m"^ in low S/N cases, no conclusions can 
be drawn with respect to the polarization angle at A -Q. 

the weight function by a sum of 6 functions. That enables us to 
discretize equation Q: 



{F*R}(4>k)^ 



2,= 1 



(3) 



where f ,■ - PiA^) for brevity, and w,- is the weight of a data 
point. We refer to the companion paper for a formal derivation 
of these equations. The RMTF is given by 



2/=i 



(4) 



One is free to choose Aq. We chose the wavelength cor- 
responding to the weighted average A^, which minimizes the 
component of the transfer function orthogonal to the actual po- 
larization vector at /Iq. 

Fig. 12 shows the RMTF corresponding to the sampling of 
A^ space that was obtained in this work. It shows the amplitude, 
the real part, and the imaginary part of the RMTF. The real part 
is the response parallel to the polarization vector at Aq and the 
imaginary part is the response orthogonal to it. The output of 
the RM synthesis procedure is a data cube with axes a, 6, and 
((>. It has complex values P - Q + iU. 

4.2. Tlie RM-cube 

A total of 126 complex polarization images were used in the 
construction of the RM-cube. A range of Faraday depths from 
-300 to +300 rad nT^ was synthesized. Beyond this range no 
features were detected. The recognition of real astronomical 
structures as well as instrumental artefacts is best done by an- 
imated scanning through the cube. Such a movie will be made 
available on our web page to help the reader in the recogni- 
tion of the observed structures. A selection of frames from the 
RM-cube is shown in Fig.0 
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RM-synthesis adds polarization images over a considerable 
range of frequencies ranging from 315-370 MHz. These im- 
ages would normally have angular resolutions that differ by a 
factor 1 .2, or a factor 1 .44 in beam area. In generating the half 
resolution image cube we therefore tapered to 1/e at a base- 
line of 1355 wavelengths, providing an almost frequency in- 
dependent beam. In principle we could also correct the im- 
ages for the different primary beam attenuations. However, 
this would cause a frequency dependent noise level which de- 
pends on radial distance. A full discussion of the various issues 
in dealing with primary beam attenuation, and the related is- 
sue of the intrinsic source em ission spectrum, is presented in 
iBrentiens & de Bruvnl ( l2005l) . 

We have also made RM-cubes with full angular resolution, 
but we will only discuss in detail the half resolution images 
which have better S/N for the generally extended emission. We 
have verified that this angular smoothing, which locally could 
lead to some beam depolarization, did not affect our conclu- 
sions in any significant way. 

The addition of 126 complex images, each with a noise 
level of about 0.7-0.8 mJy, resulted in an RM-cube with rms- 
noise levels going down to about 70 fiiy beam ' in Stokes Q 
and U. However, these levels are only reached at the edge of the 
field and/or at high RM values (> 100 rad m"^) where instru- 
mental problems have been 'wound up' sufficiently to be left 
only with the (uncorrected) noise in the individual channels. 

5. Description of observed structures 

Fig. shows several representative frames from the total po- 
larization RM-cube. Most of them are separated by about 
9 rad nT^, which should be compared to the resolution in RM- 
space of 15 rad m"^. These images have not been corrected for 
the reduced off-axis sensitivity in order to preserve a uniform 
noise level across the images. As discussed above, the images 
produced using the RM-synthesis technique utilize an effective 
bandwidth of about 50 MHz yielding a sensitivity of better than 
100 juJy beam \ This is more than an order of magnitude be- 
low the total intensity noise level and with a peak brightness 
in Stokes / of 22 Jy (in the half resolution image) this repre- 
sents a formal dynamic range of about 200 000: 1 . It is clear, 
however, that the dynamic range varies enormously in the RM- 
cube. Thus in addition to the real (astronomical) signals from 
cluster synchrotron emission, the polarized images as well as 
the RM-cube images show many structures and patterns that 
have an instrumental origin. We will therefore begin with a 
description and, where known, an explanation of these instru- 
mental artefacts. This will be followed by a description of the 
astronomical features. 

5.1. Instrumental artefacts 

We have identified several prominent patterns in the polariza- 
tion and RM-cube images that clearly have an instrumental ori- 
gin. In recognizing instrumental artefacts it is important to re- 
alize that they are a combination of uv-plane and image-plane 
effects. Multiplicative errors in the uv-plane result in a convo- 
lution with an error pattern in the image plane. The strongest 
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Fig. 3. Faraday spectrum at the location of the east (E) and 
west (W) lobe of the extended radio source at a ^ 3^10", 6 ~ 
42°50', revealing the instrumental polarization resonances for 
strong off-axis sources. 



effects are therefore associated with the strongest sources, al- 
though they can extend over a substantial part of the image 
plane. The off-axis effects, however, also have a multiplicative 
character resulting in a spatially dependent convolution with 
an error pattern. Again the strongest sources show the largest 
effects. 

Radio emission observed away from the telescope pointing 
direction of reflector antennas is generally polarized. The com- 
ponent that is uniform across the field of view is small (a few 
%) and is calibrated away using standard procedures (Section 
\3.2l . The instrumental polarization (in Stokes parameters Q, 
U and V) in parabolic dishes, however, rap idly increases with 
distance from the optical axis ('NaDiei'1999'). In the case of the 
WSRT it increases to about 20% at the -10 dB power point 
with a typical clover-leaf pattern. Stokes Q is positive in the 
E-W direction and negative in the N-S direction. The Stokes U 
pattern is rotated counterclockwise relative to Stokes Q by 45°. 

The off-axis polarization causes spurious signals at the lo- 
cation of the many strong sources in the field (see Fig.0}. The 
inevitable small spread in off-axis polarization levels for differ- 
ent telescopes in the array also causes weak ring-Uke patterns 
(see Fig.|4}. Furthermore, the small variations between differ- 
ent 12 h observations also cause weak residuals at multiples of 
the grating lobes radius (about 40', and frequency dependent). 
Intersection of these grating lobes then gives rise to a number 
of weak spurious features. 

At RM=0, all contributions add up coherently, thus if the 
off-axis polarization would be independent of frequency the re- 
sponse in the RM-cube should rapidly diminish at large RM 
values. The observed response as a function of RM indeed 
drops fairly rapidly but not as fast as expected. This is due to 
the fact that the off-axis polarization of the WSRT, in Stokes 
U and to a lesser extent Stokes Q, also reveals a very strong 
frequency dependence with a period of about 17 MHz.^ 



- This instrumental attribute has been known to exist for a long time 
but its origin is still not fully understood. It is believed that there is a 
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Right Ascension (J2000) Right Ascension (J2000) 

Fig. 4. Derotated Stokes QdX?? - 0.759 for the foreground Fig. 5. Polarized intensity ||P|| at 0= 52 rad m"^. The grey scale 

emission at a Faraday depth (p - 6 rad mT^. The intensity scale covers a range from 0.3-2 mJy beam ' rmtf"' . 
saturates at +1.5 mJy/beam (white) and -1.5 mJy/beam (black). 



The 17 MHz period causes peaks in the RM 'spectrum' 
at values of about +42 rad mT^ and -42 rad m"^, in addi- 
tion to the emission at rad voT^. These 'resonances' can be 
clearly discerned in intense sources located far ofF-axis when 
scanning through the RM movie. The elongated source in the 
north-western upper corner of the image (Fig. ^ shows this 
very clearly (Fig.O. Detailed scrutiny of the Faraday spectrum 
over the face of the source suggests that the part of the emission 
around is +5 rad m"^ may in fact be partly due to intrinsic po- 
larization from the source. The predicted Galactic value at this 
location (see Section^] Fig. II U is about -10 rad voT^. 

A third instrumental artefact in the images is a faint pattern 
of radial stripes emanating from 3C 84, by far the strongest 
radio source in the field. The intensity of these spikes rapidly 
dies out when we move away from position angles of ±90°, 
giving them the appearance of 'whiskers'. Their peak intensity 
is about 0.5 mJy beam"' but they are more typically present 
at 0.1-0.2 mJy beam These whiskers are best seen in Fig.|5l 
but upon closer inspection they can also be seen in several of 
the frames of Fig.0(e.g. RM= -24, 21 and 51). Their origin is 
still being investigated. They were present only in the original 
Stokes U (and V) images but the rotation of the polarization 
vector, after de-applying the ionospheric Faraday rotation, as 
well as the rotation inherent to the RM synthesis process itself, 
leads to their appearance in every RM frame. 

Finally, we also detect polarized grating lobes from the 
source Cas A in the input Q and U image cubes. However, due 



component due to a standing wave pattern between the focus box and 
the dish (separated by 8.75 m) and a component due to scattering off 
the four legged support structure of the focus box. 



to the great angular distance of this source from the field centre 
(» 40°) its chromatic grating response moves rapidly across the 
field of view as we move in frequency. In the final RM-cube, 
these instrumental features therefore decrease to a level well 
below the noise. 

5.2. Astronomical signals 

Once the instrumental features are recognized, the polariza- 
tion images and the RM-cube frames are seen to be filled with 
highly significant signals of clearly celestial origin. They cover 
a wide range of Faraday depths (p between and 90 rad m~^. 
Broadly speaking there seem to be two types of emission: 

1 . diffuse structures with slowly changing polarization angles 
on scales of the order of several tens of arc minutes: these 
have < < 15 rad nT^ 

2. distinct large structures having sizes of the order of a degree 
and granularity in the polarization angle on scales of the 
order of a few arc minutes: these have 30 < < 90 rad m"^ 

There appears to be no significant emission between 15 and 
30 rad nT^. 

The first type of diffuse emission can be seen best in the 
third frame of Fig. This diffuse emission appears to be spread 
over an area significantly wider than the primary beam (2^2- 
2°.6) although it does decrease in brightness, as any astronom- 
ical emission ought to. We emphasize this because the high 
Faraday depth emission is notably more confined to the central 
area. Fig. |3 shows de-rotated Stokes 2 at a Faraday depth of 
+6 rad m"^. It illustrates the large coherence of the polariza- 
tion angle across these patches. 
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Fig. 6. Faraday depth spectra of two regions of Galactic fore- 
ground emission. The upper spectrum shows a region slightly 
west of the bar. The lower spectrum shows a foreground patch 
near the lens. 



6. Galactic foreground emission at low (p 

The large-scale diffuse polarized emission at low <p is very 
similar to the features seen in previous WSRT studies at 
92 cm. This component has been shown to be produced in 
the relatively local Galactic mediurn, probably well within 1 
kpc from the Sun jHaverkorn et al.ll2003at rHaverkornll2002t 
IWieringa et al.lll993 ). This diffuse and complex polarization 
structure is believed to result from the line-of-sight superpo- 
sition of intrinsically highly polarized emission and spatially 
as well as depth dependent Faraday rotation by the interstel- 
lar (Galactic) magneto-ionic medium. Haverkor n et al.l ll2003bl) 
found that the rotation measure of this Galactic foreground in 
Auriga (/ = 161°, b - +16°) is typically between -17 and 
H-10 rad nr^, with an average of -3.4 rad m"^. The values we 
find towards the Perseus cluster at similar / (I « 150°) but oppo- 
site b (b -13°), are very similar. They cover the range from 
Oto 12 rad m 2. 

Fig. |6l displays a few typical Faraday depth 'spectra' of 
Galactic foreground patches, in areas where there is no signifi- 
cant emission due to (side lobes from) instrumentally polarized 
emission. It is also evident that at least part of the foreground 
is 'extended' or 'thick' in Faraday depth. Both spectra suggest 
that part of the positive contribution to the Faraday depth oc- 
curs inside the medium that produces the foreground emission. 
One should nevertheless not be misguided by the height of the 
peaks in a Faraday spectrum. The various artefacts in the indi- 
vidual Q and U images that were described above occasionally 
conspire to produce spurious peaks and one should therefore 
inspect the individual Q and U images at the location of a peak 
in a RM-movie to verify the reality of the RM-feature. Spatial 
continuity is a powerful aid in ascertaining the reality of most 
features, and this is best done in the RM-movie. 

We find no exceptionally bright foreground features di- 
rectly in front of 3C 84 that could be associated with the low 
surface brightness polarized emission that was detected in very 



high dynamic range 21 cm WSRT observations (^Siibr inJl993t 
de Bruvn 1995) . We will return to this in Section lOl 

7. Background structures at high 

7.1. Description of structure in the RM-cube 

The second type of diffuse polarized emission is much richer 
in spatial structures. Note that we will only discuss the polar- 
ized intensity and defer to a future paper the discussion and 
analysis of the rich structure in the polarization angle distri- 
bution. The description starts with the fifth frame of Fig. 
{(f> = 30 rad m"^) which shows a weak, front-like structure run- 
ning from a ^ 3H6"',6 ~ 40°24' to a ^ 3H0'",6 ~ 41°36'. 
The next frame {<p - 42 rad m"^) displays a stronger linear 
feature with a slightly different position angle. A bright cir- 
cular 'doughnut shaped' structure, with a diameter of about 
7' develops at position a ^ 3*'15'"35%5 ~ 41°42:3. Fig. H 
shows Faraday depth spectra of two lines-of-sight through the 
doughnut. The dashed spectrum shows the high significance 
of the doughnut. The peak at negative Faraday depth in the 
solid spectrum is caused by a whisker. The bright peak near 
<p = +10 rad m - is due to the Galactic foreground. 

Frame number seven - 51 rad m"^) shows a spectacular 
lenticular feature southwest of the doughnut. It is even better 
seen in Fig.|5l which is an enlarged frame from the RM-cube 
at <p - +52 rad m"^. The position angle of the lens is very 
similar to the position angle of the linear structure in the pre- 
vious frame. The lens is roughly 40' x 20' in size. At the dis- 
tance of the Perseus cluster, it translates to 1 x 0.5 Mpc. Two 
Faraday depth spectra through the lens are shown in Fig. |9l 
The two peaks at negative (p in the solid spectrum are caused 
by instrumental whiskers. The emission from the lens peaks at 
4> ~ 52 rad mT'. This frame and the next one, frame 8 at = 60 
rad mT^, also show bright extended emission with a polarized 
intensity of about 1 mJy beam"' approximately 40' north of 
3C 84. This patch is roughly 15' x25' and is located north-west 
of the steep spectrum tail of NGC 1265 dSiibring & de Bru'viil 
Il998l) . 

The lens and doughnut slowly fade away in the frames that 
follow and polarized emission disappears from the western part 
of the field. Frame nine shows mottled emission across an area 
of 2° X 1?5 centred around the area between NGC 1275 (3C 84) 
and NGC 1265. The tenth frame (0 = 78 rad m -), is the last 
frame that shows significant emission. The horizontal bar al- 
most 1° due north of 3C 84, at a 3''20™, 6 ^ 42°25', is very 
conspicuous at these values of <p. Two Faraday depth spectra are 
shown in Fig. [TO] The emission fades away towards the cluster 
centre. Beyond (p ^ \Q0 rad m"^ no structures are detectable at 
2'-3' resolution. 



7.2. Where is the total intensity counterpart I? 

We have not yet detected any of the extended polarized fea- 
tures in our RM-cube in the total intensity image. This was 
rather surprising because we have no doubt that this emission 
is due to the synchrotron process. There may, however, be a 
rather simple explanation. The sensitivity in the Stokes / im- 
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Fig. 7. Characteristic frames of polarized intensity (\\P\\) from the RM-cube. The Faraday depth of each frame is specified in the 
top left comer in rad m"^. 
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Fig. 8. Faraday depth spectra of two lines-of-sight through the Fig. 10. Faraday depth spectra of two lines-of-sight through the 
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straight bar north east of NGC 1265. 
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Fig. 9. Faraday depth spectra of two lines-of-sight through the 
lens. 



age is about a factor 10-15 poorer than in the polarization im- 
age: 1.5 mJy beam"' as against 0.1 mJy beam"' (both numbers 
now refer to the 2' x 3' resolution image. However, the bright- 
est features in the polarization images, e.g. the iens' and the 
'doughnut', have a peak brightness of about 1 mJy beam"' . So 
the lack of detection in / implies, at face value, a polarization 
percentage of at least 50%. Although this is not impossible for 
synchrotron emission, these percentages are getting 'uncom- 
fortably' high. Seve ral 'relic' structures in clusters of galaxies 
(lEnfilin etalJI 19981) have polarization percentages of 20-30% 
although these all refer to shorter wavelengths where depolar- 
ization is less of an issue. 

Apparent polarization percentages well above 100% are 
quite normal in t he case of Galactic foreground polariza- 
1993: 



tion (iW ieringa e t al. 1993: Haverkorn et al. 2003a; Haverkorni 
12002! lOaensler et alJl200lh . Thev are due to the fact that the 
Galactic foreground Stokes I is very smooth and largely re- 
solved out due to the lack of sufficiently short interferometer 
spacings. The emission in Stokes Q and U, on the other hand, 
has high spatial frequencies because they can be both positive 
and negative (cf. Fig.|3. Fig. ^2 to which we will return later, 



shows that the distribution of polarized intensity indeed has a 
significant large-scale component. If the total and polarized in- 
tensity would be distributed similarly, the shortest spacings of 
the WSRT will have only partly sampled the signals. This was 
exacerbated by the fact that in two out of six usable frequency 
bands one of the crucial telescopes (RT9 and RTA) providing 
the shortest 36 m spacing were not functioning properly. The 
true polarization percentage of the polarization structures could 
then easily be a factor 2-3 lower which would decrease the ap- 
parent percentage polarization from >50% to >20%. The con- 
clusion, however, remains that the signals must be highly po- 
larized. We will return to this aspect below. 

8. Discussion 

8.1. Location of high Faraday depth features 

There are two possible locations for the source of the emission 
at (p > 30 rad m"-. The first possibility is that here too we 
are dealing with Galactic synchrotron emission. Potential sites 
are the region between us, and possibly including, the Perseus 
arm which is at a distance of 2-3 kpc. The second possible 
location is the Perseus cluster of galaxies, which was the central 
target of these observations. We favour the second option for a 
number of reasons: 

- the high structures seem to be largely confined to and 
approximately centred at the Perseus cluster, contrary to the 
low (p emission; 

- there is a systematic decrease in the small scale structure 
in Stokes Q and £/ as a function of distance to the cluster 
centre. 

Initially we believed there was a third argument for an associ- 
ation with the Perseus cluster: when scanning the RM-movie 
there is a very suggestive systematic convergence of the higher 
RM values towards the centre of the Perseus cluster. However, 
we could not exclude that this was caused in part by an anomaly 
in the Galactic foreground RM distribution. 
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Ijohnston-Hollitt et alJ (l2004l) have constructed an interpo- 
lated map of the RM distribution of a sample of over 800 ra- 
dio sources, from which we may estimate the Galactic RM 
contribution in any direction. The expected RM near the 
Perseus cluster is somewhere between H-10 and -1-20 rad nT^. 
The RM expected from the large-scale Galactic field should 
be n egative in this area. The positive values inferred by 
Ijohns ton-Hollitt et al. (2004) may therefore indicate an RM 
anomaly south of the Galactic Perseus arm. 



8.2. Evidence for a Galactic RM gradient 

To investigate this further and put our suspicions about fore- 
ground contamination to the test we proposed a new series 
of observations in August 2004, when most of this work had 
already been finished and written up. Using the WSRT at 
21 cm we conducted a series of snapshot observations of 15 
polarized (background) radio sources in the direction of the 
Perseus duster The sources were selected through the NVSS 
( ICondonetalll998l) and should have at least a few mJy of po- 
larized emission. Using the 160 MHz wide WSRT band this 
appeared to be sufficient to allow a RM determination accu- 
rate to better than a few rad m"^ for all 15 sources after 30 m 
of integration per source. After excluding three sources with a 
complicated brightness morphology and multiple, widely dif- 
ferent RMs, the remaining 12 sources showed a smooth gradi- 
ent across the field of view. The results are shown in Fig. II Hand 
include a radia l basis-function interpolation of a smooth fore- 
ground screen jCarr et alJ200ll) . Because the uncertainty in the 
RM determination of each source is at most a few rad m"^, the 
gradient is very significant. We also infer that the scatter due 
to a RM contribution from the medium inside, or surrounding, 
the radio sources can be at most 5-10 rad m"^; this is consistent 
with previous studies of RM diff'erences in double radio source s 
jSimonetti etalJl984HSimonetti & Cordesll986HT,eahvll987h. 

This smooth Galactic foreground RM gradient should be 
compared with the observed spatial distribution in the diffuse 
high Faraday depth emission features. To do this we proceeded 
as follows but we begin with noting that there is no such thing 
as the RM or the Faraday depth of a pixel: indeed there can 
be multiple regions of emission at different Faraday depths 
along the line-of-sight (for an extensive discussion of such 
situat ions, which are quite normal at low frequencies, we re- 
fer to lBrentiens & de Bruviil ll2005l) V Within the RM-cube we 
masked out the emission at 3 < (p < 19 rad m"^ which is 
very clearly dominated by the Galactic foreground. At every 
pixel in the RM-cube we then determined at which Faraday 
depth the polarized intensity \P\ was maximal. This range in- 
cluded ^ - Q rad m"^, the value where all frequency indepen- 
dent instrumental problems accumulate. This guaranteed that 
the emission that would be assigned a Faraday depth other 
than zero, cannot be caused by RMTF side lobes of instru- 
mental problems. The regions where the latter is manifestly 
the case are the sources at the edges of the field. The result 
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Fig. 11. Faraday depth of the maximum |P| along the line-of- 
sight, excluding range 3-19 inclusive. Pixels with a peak po- 
larized flux density less than 0.7 mJy per beam have been 
masked out. Ove rplotted is a radial basis function interpolation 
jCarr et alJ200ll) of the "foreground RM" determined using the 
12 sources indicated with filled circles. The thick line has a RM 
- rad m"^, and contours are given in steps of 20 rad nT^ with 
positive values increasing to the east. 

of this exercise is shown in Fig. II If . We imposed a cut at 
|P| < 0.7 mJy beam"' before including pixels. 

Fig. ^Jrsveals that across the 2° diameter area with high 
(p emission there is a clear Galactic RM gradient of about 
60 rad m"^, which is approximately equal to the observed range 
of (f) - 30-90 rad nT^. In other words, the high (/> emission 
shows a fairly uniform positive excess of about 40-50 rad m"^ 
relative to the interpolated Galactic RM. For all but two of the 
background sources the line-of-sight avoids the Perseus clus- 
ter, which we define as the area out to which significant X-ray 
emission can be seen, a radius of about 1°. We may therefore 
assume that the RM of these background radio sources rep- 
resents the line-of sight-integrated RM of our Galactic fore- 
ground magneto-ionic medium. This result holds the key to the 
interpretation of the high </> emission as we will now proceed to 
argue. 

For the sake of the argument let us concentrate on the fea- 
tures on the western side of the cluster near what we have called 
the shock, the lens and the doughnut. The diffuse organized 
Galactic polarized emission has a RM of about -I-IO rad m"^. In 
order for the total, integrated, RM to the background sources to 
end up at about rad nT^ there must be another small negative 
contribution of about -10 rad m"^ from somewhere along the 
line-of-sight in our Galaxy. This seems not implausible. The 
RM-cube around ~ rad m"^ is rather messy due to the 
remaining instrumental polarization from the extremely bright 

^ The on-line version of this article contains a colour version of this 
figure. 
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and complex cluster radio sources. Now let us suppose that the 
'screen' responsible for the 40-50 rad m"^ contribution is also 
due to our Galaxy and extends across the full field of view 
shown in Fig.^2(the sudden termination at a distance of about 
1°5 from the pointing centre must then be attributed largely to 
the primary beam attenuation, which is unlikely but not impos- 
sible). We would then require another screen, but now with a 
Faraday depth of (p - -40 rad mT^ to compensate for the (p - 
+40 rad vnT^ screen. We see no evidence in our RM-cube for 
any emission from such a screen. 

Even stronger arguments can be brought forward for the 
case that the (p - +40 rad mT^ emission would result from a 
discrete cloud in the Perseus arm of our Galaxy that just hap- 
pened to cover the Perseus cluster Assuming a magnetic field 
of 1 juGauss, a depth of about 60 parsec (equal to a lateral scale 
of V.5 at 2.5 kpc), a RM of 40 rad mT^ requires a (smooth) 
electron density of 0.8 cm"^. Assuming roughly 20% Ha ex- 
tinction, this density and path length would generate an Ha 
surface brightness of about 13 Rayleighs. The total Ha surface 
brightness towards t he Perseus cluster in the WHAM survey 
llHaffneret alJ 1200 3) is about 4 Rayleighs. Furthermore, the 
integrated Ha emission is rather smooth towards the Perseus 
cluster, ranging in surface brightness from 3-5 Rayleighs over 
an area of 3° x 3°centred on the Perseus cluster There is no 
excess at the location of the cluster The discrete Perseus arm 
model for the origin of the excess +40 rad mT^ polarized emis- 
sion therefore appears to be untenable as well. 

A similar argument can be given for the emission to the 
north and east of the cluster centre. The emission in this region 
shows a significant drop in surface brightness well before the 
primary beam attenuation sets in. This would have to be a pe- 
culiar coincidence in the case of a Galactic foreground origin. 

When we add the above arguments to those given previ- 
ously, the spatial coincidence with the cluster and the granular- 
ity of the polarization angle structure, we are led to the conclu- 
sion that the high (p emission indeed must be associated with 
the Perseus cluster of galaxies. 

8.3. Intracluster or peripheral emission? 

Having argued that the high ^ emission is associated with the 
Perseus cluster we may then ask the question: where in the clus- 
ter does the emission originate? Does it originate from within 
the cluster or is it emitted from the periphery of the cluster? At 
this point we should also realize that this question really has 
two aspects: 

1 . Where does the emission come from? 

2. Where does the Faraday rotation occur? 

Fig-Elis a sketch that is helpful for the following discus- 
sion. It shows where we believe the emission and rotation oc- 
cur in the Perseus cluster From the fact that most of the po- 
larized structures are unresolved by the RMTF it appears that 
the m edium is not Faraday-thick (see Brentiens & de Bruvn 
( l2005b V This means that is the largest fraction of the Faraday 
rotation of 40 rad mT^ occurs between the emitting region and 
us. 



rotation 




Perseus cluster 



> Polarized emission 




Faraday rotation (+40 rad m ^ ) 
I Earth 

Fig. 12. Sketch of our proposal for the situation in the Perseus 
cluster A top view of the distribution in the Perseus cluster of 
various emitting areas and Faraday rotating areas. 



At this point it is relevant to recall that there are two ex- 
tended and polarized radio sources within the Perseus clus- 
ter, NGC 1265 and IC 310, in the area where we find diffuse 
polarized emission. Their head-tail morphology indicates that 
they are most likely located within the denser gaseous part 
of the Perseus cluster (Miley et al. 19 721). The RM across the 
brigh t part of the tail of NGC 1265 lo'Dea & Owenlll986l 
Il987l) shows a scatter of 20 rad m ^ around a mean value 
of 45-50 rad mT^. The scatter is probably due to the magne- 
tized plasma in the cocoon surrounding and mixed within the 
twin tails of this head-tail source. The 'expected' value due to 
our Galactic foreground at the location of NGC 1265 is about 
25 rad m"^ (cf. Fig. II 1> . The difference of about 20 rad m"^ 
must then be attributed to the medium in the line-of-sight be- 
tween the head-tail source and the edge of our Galaxy: i.e. the 
peripheral region of the near side of the cluster. 

We also detect faint but highly significant polarized emis- 
sion from a region close to the 'head' of IC 310 at a value of 
(p - +80 rad m"^ (Fig. I16> . This polarized emission probably 
originates from a region about 1' downstream from the head of 
IC 310 where significant polarized emissio n was also detected 
at 610 MHz bv lSiibring & de Bruvnl(ll998h . This is sufliciently 
far away from the galaxy that we may assume that the Faraday 
rotation is not due to the ISM associated with IC 310 itself. 
The expected Galactic RM value (cf. Fig. II 1> at this location is 
about +20 rad mT^. The difference of +60 rad mT^ could imply 
that IC 310 is located deep within the cluster. 

Note that the excess Faraday depth shown by NGC 1265 
and IC 310 is consistent with the average cluster RM excess of 
less than 50 rad mT^ found by Clarke et al. (2001) at the pro- 
jected distance of 700 kpc of these two head-tail radio sources 
from the centre of the Perseus cluster The polarized emis- 
sion having (p- 40 rad voT^ is probably located on the periph- 
ery on the near side of the cluster at about the same depth as 
NGC 1265. IC 310 may well be located behind the diffuse po- 
larized emission. 
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Fig. 13. Thomson scattering geometry. 



8.4. Thomson scattering within the Perseus cluster? 

As discussed in the introduction the initial goal of the observa- 
tions was to characterize the possible contamination of Galactic 
foreground pola rization to the observed 21 cm polarization 
jde Bruvnlll995h which was observed to straddle 3C 84 on an 
angular scale of 30' and believed to be due to Thomson scat- 
tering. In a future paper we plan to present the original 21 cm 
observations augmented with a wider field study. 

Could the polarized emission that we detected at wave- 
lengths of 81-95 cm be indeed Thomson scattered radia- 
tion from 3C 84 or its direct surroundings (ISvunvaevlll982t 
[wise & Sarazin 1990, 1992)^? 

In order to quantitatively answer that question, we have 
made a simple model of the Perseus cluster. The redshift of 
the Perseus cluster is z = 0.0167 (Struble & Rood 1999). We 
assume i/o = 72 + 2 km s"' Mpc"' ( S oereel et al.l2003ll) which 
gives a distance to the cluster of 69.5+2 Mpc. We have adopted 
a spherically symmetric electron distributio n. We used the 
deproj ected radial dependence determined bv IChurazov et alJ 
1I2OO3I): 



3.9 X 10 



-2 



4.05 X 10- 



V80 kpc/ ^ ^ \ 280 kpc) 



0.87 



(5) 



In order to estimate the radio radiation, we summed the 
contributions of the 3C 84 30" component and the halo. We 
ignore the very bright core, because it has an uncertain his- 
tory and its high flux density is largely due to an outburst that 
only started in the late 1950's jO'Dea et al.ll984l) . Its Thomson 
echo cannot therefore have spread by more than a few tens of 
parsecs. The values for t he lumiriosity fr om the extended com- 
ponents are taken from ISiibrind (Il993l) after scaUng them to 
our adopted distance. Both source components were situated in 
the centre of gravity of the electron population. We modelled 



Please note that lWise & Sarazir] il990l) made a small error in their 
equation (1), which m akes it incompatible wit h our equations J10> and 
irn . Equation (1) of IWise & SarazirJ <1990l) must be multiplied by 
l/4;r in order to correct this. Th e same error is also present in equation 
(2.1) o fl Wise & Sarazi hence their estimates are one order of 

magnitude too optimistic. 
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Fig. 14. Computed Thomson scattering halo surface brightness 
at 92 cm based on the radio flux density of the 3C 84 halo and 
30" component. Total intensity. 
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Fig. 15. Computed Thomson scattering halo surface brightness 
at 92 cm based on the radio flux density of the 3C 84 halo and 
30" component. Polarization fraction. 



the source as a single unpolarized point source and assumed an 
isotropic radiation field. 



P(v') = 10.2 X 10- 



24 



9.83 X 10 



24 



333 MHz 



333 MHz 



(6) 



-1.2 



WHz 



where v' = (1 -H z)v is the rest frame frequency of the Perseus 
cluster. 

The geometry of the computation is given in Fig. The 
radio source is located at the origin (O). The emission is scat- 
tered at location (b, x). The total intensity of the scattered emis- 
sion is given by the line integral 



f^- P((l+Z)v) dcr 
/(v) = J_^ ^^^ne(r(.,^))dx, 



(7) 
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where 



dQ 
and 



(8) 



(9) 



llSvun vaev'l982':'R vbicki & L ightman' 1979V We can write 
equation (0 in terms of b and x as depicted in Fig. 1131 



/(v) = K(v) J n^{r)— -^d;*: W Hz"' mT^ Sr"', (10) 



where b is the impact parameter, x is the path length measured 
from the midplane of the electron distribution, r - y/b^ + x?- 
is the distance from the radio source, L = V^max ~ 

is the 

integration boundary, and 



3f((l +z)v)o-T 
64;r2(l 



W m-" Hz ' Sr" 



(11) 



independent of position. We used SI units throughout this com- 
putation. Assuming there is no Faraday rotation in the clus- 
ter medium, which of course is unrealistic at a wavelength of 
92 cm, we may also compute the maximum degree of polariza- 
tion by 



/p(v) = K{v) 



(b^ + x2)- 



-d.xWHz"' m"^ Sr"'. (12) 



The results are shown in Fig. 1141 and Fig.^] Fig.ll4lindi- 
cates that at small projected distances (< 30') from the cluster 
core it is, in principle, possible to detect the Thomson scatter- 
ing emission even at 92 cm wavelength. For example, at 25' 
distance from the cluster centre, the expected total intensity is 
approximately 50 /iJy 2' x 3'beam"'. Assuming 50% polariza- 
tion, this amounts to a polarized intensity of 25 yuJy beam"'. 
Barring dynamic range issues, it will however be very diffi- 
cult to separate this scattered emission from true cluster halo 
emission in Stokes /. The much deeper sensitivity reached in 
polarized intensity (70 yuJy beam ') affords a better chance if 
we can correctly estimate the complications from Faraday ro- 
tation within the scattering medium. So although we can not 
exclude that part of the mottled polarized emission close to 
3C 84, at the location where we previously observed significant 
polarized emission at 21 cm, is due to Thomson scattering, we 
will not discuss it here any further and will concentrate on the 
brighter polarized features. 

We believe that neither the large structures that we detected 
in the western part of the field, nor the bar feature to the north, 
can be explained by Thomson scattering. The high percentage 
of polarization inferred for the diffuse structures are consistent 
with Thomson scattering (without significant de-polarization), 
but they appear far too bright for their projected distance from 
3C 84. For example, to explain the observed polarized surface 
brightness of the lens (approximately 1 mJy per 2' x 3' beam) 
by Thomson scattering, either the source luminosity 3 to 6 mil- 
lion years ago should have been a factor of 100 to 500 higher. 
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Fig. 16. Faraday depth spectrum near the 'head' of IC 310. A 
highly significant feature at Faraday depth = 80 rad nT^ can 
be detected. The peak around (p - rad nT^ is due to the in- 
strumental polarization of IC 310, which has a peak in total 
intensity of more than 500 mJy. 
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Fig. 17. Faraday depth spectrum a line-of-sight through the 
curved blob due north of NGC 1265. 



or the electron density at 1 to 2 Mpc distance from the cluster 
centre should be at least 5 x 10"^ cm""*. Both options, or a com- 
bination of a higher luminosity and a higher density, appear to 
be implausible. However, we note that the current core lumi- 
nosity is 5x higher than the luminosity assumed in the com- 
putations and the uncertainty about the long term average of 
the core emission will remain one of the uncertain issues with 
the Thomson scattering interpretation. In a future paper we will 
return to the issue of the separation of Thomson scattered emis- 
sion from other Perseus cluster contributions using both exist- 
ing and new 21 cm, 92 cm, and 200 cm WSRT observations. 

8.5. Structure formation shocks and AGN bubbles 

We will now briefly discuss the implications of what we 
consider to be the most likely explanation for the nature 
of the diffuse high Faraday depth structures centred on the 
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Perseus cluster. A more in-depth study will b e presented 
in future papers. We will fol low lEnBlin et al.l {l998) and 
lEnBlin & Gopal-Krishnal (l200lh and interpret the emission as 
resulting from (re-energized) relics and emission associated 
with shocks in the large-scale structure formation gas flow. 

Cosmological simulations of structure formation have im- 
proved dramatically in recent years. Nearly all simulations 
indicate that structure in the Universe formed hierarchically. 
Small clumps collapsed, larger ones started accreting smaller 
ones and grew even larger A consequence of this type of 
structure formation is the presence of large-scale gas flows, 
which fo rm huge, mul t i megaparsec scale shocks at their inter - 
sections teurns 'iggs'; 'Ouil is et alJll998l: iMiniati et alJbOOOl) . 
lEnBl in et al. ( 1998) discuss the possibility that cluster relic 
sources, large highly polarized radio sources in the outskirts 
of galaxy clusters, could trace interface shocks between clus- 
ters and the super cluster filaments accreting onto them. The 
electrons in flie relic would be re-energized b y diffusive shock 
acceleration. 'EnBl in & Gopal-Krishnal ( l200ll) developed an an- 
alytic formalism for an alternative process. They propose that 
cluster rehcs are bubbles of magnetized plasma, released into 
the IGM by AGN. The electrons in such a bubble loose en- 
ergy through adiabatic expansion of the bubble. Once fully de- 
tached, the bubble enters a 'hibernation' phase in which it can 
stay for several Gyr; the only energy losses suffered by the rel- 
ativistic electrons are those due to Inverse Compton losses off 
the 2.7 K background radiation. T he bubble is now buoyant in 
the IGM, due to its lower density ('Bruggen'2003'). It is finally 
revived when the bubble is adiabatically comp ressed by a l arge- 
scale structure formation shock tenBhn & Bruggenl2002l) . 

The morphology of the long linear structures that appear in 
frames 5 and 6 of Fig.0are indeed reminiscent of shock fronts. 
This structure would then have to be located at the periphery 
of the cluster As argued in section lOl it would probably have 
to be located on the near side of the cluster to prevent signifi- 
cant beam depolarization following the passage of the radiation 
through the Perseus cluster. The good alignment of the front 
(frame 6) and the lenticular structure (frame 7) argues that they 
are co-located. Do we perhaps observe the actual compression 
of a huge (1 x 0.5 Mpc) hibernating bubble by a large-scale 
structure formation shock? 

The Perseus cluster is situated at the eastern end point of the 
Perseus-Pisces super cluster The fronts and the lens are in the 
area where one would expect large-scale structure formation 
shocks from general structure formation simulations. The fact 
that the front is so straight over a distance of more than 3 Mpc 
is not impossible; many of the fronts in these simulations show 
that behaviour EnBlin et al. (2001) report a tentative detection 
of a cosmological shock wave that overruns one of the radio 
lobes of NGC 315. NGC 315 is positioned at the intersection 
of several sub filaments of the Perseus-Pisces super cluster 

The structure of the doughnut strongly resembles the struc- 
ture and topology of the magnetic field in a pre-shock bub- 
ble, as simulated by .EnBlin & Briiggen (.2002 ). Assuming the 
doughnut to be situated at the distance of the Perseus cluster, it 
measures 250 kpc across, consistent with the hypothesis of an 
AGN bubble. 



The blob that is visible in frame 7 of Fig. north of 
NGC 1265, is also a good bubble candidate. It has a slightly 
curved shape. Its angular size is approximately 30' x 15', which 
corresponds to 750 x 375 kpc at the distance of the cluster It 
is striking that it is located just north of the steep spectrum tail 
of NGC 1265 and has roughly the same curvature. As can be 
seen in Fig. [n] its Faraday depth is approximately 50-60 rad 
m This is close to the value observed in the main part of 
the tail by O'Dea & Owen ( 1986). Could this structure be a de- 
tached bubble associated with a previous phase of activity of 
NGC 1265? The projected separation of about 15' translates to 
about 300 kpc or 300 million years for a transverse velocity of 
1000 km sec We have also searched for pola rized emission 
in the steep spectrum tail of NGC 1265 ( Siibri n^ & de Bnivi] 
Il998h . but did not find any significant (^ 1 mJy beam ') peaks. 

The bar structure still further north is rather straight, but 
is both narrower and shorter than the two fronts on the west 
side of the cluster . Such structures have been discussed by 
lEnBhn & Bruggenl(l2002) . With the current limited information 
on the detailed properties of these regions it is inappropriate to 
start making a detailed comparison. 

9. Remaining puzzles and future work 

Through the discovery of highly polarized, very extended 
and very low surface brightness emission associated with the 
Perseus cluster we believe we have uncovered a rich diagnos- 
tic tool. This tool has the potential to study the magneto-ionic 
medium, the relativistic plasma content, the magnetic field 
topology and strength and the interaction with LSS -formation 
related inflow into clusters of galaxies. iKeshet et alJ J2004i) 
predict that cosmological shock waves should be easily de- 
tectable with the Square Kilometer Array (SKA) and the Low 
Frequency Array (LOFAR) at frequencies below 500 MHz. 
They suggest that the signal could already be marginally de- 
tectable with current instruments and they are probably correct. 

In future papers, following the reduction and analysis of 
several recently acquired new full polarization datasets on the 
Perseus and other clusters, at wavelengths from 21 cm to 
250 cm, we will address some of these issues in consider- 
ably more detail. The detection of similar structures in other 
clusters would support our current hypothesis for their ori- 
gin. In order to minimize complications due to high Faraday 
depth foreground polarization confusion, these clusters are lo- 
cated at high Galactic latitude. We have chosen the Coma and 
Abell 2256 clusters. Recently, filame ntary structures have been 
observed in the Abell 2255 cluster dCiovoni et al.ll2005h . The 
surface brightness of the structures found in Abell 2255 is 
about 1 to 2 orders of magnitude higher than the structures we 
found in or near the Perseus cluster We expect that using RM- 
synthesis one could trace these structures well into the large- 
scale filaments of the cosmic web. 

On the observational side we can identify the following 
tasks that lie ahead of us: 

- Image the corresponding total intensity and determine the 
spectral properties of the emission. 
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- Determine the polarization percentage and structure over 
an even wider range of frequencies. This may allow us to 
unravel the internal plasma density and intrinsic magnetic 
field structure. 

- Determine where the boundaries of the interaction between 
cluster and large-scale structure flows in the intergalactic 
medium are located. 

- Extend the sample of clusters 

- Combine the current data with data taken at different ob- 
serving frequencies. 

The last point is necessary in order to derive the magnetic field 
direction in the plane of the sky. For this one needs the polar- 
ization angle at A - Q, xo- Unfortunately the Icr error on the 
derived Faraday depth is too large (more than 1 rad m"^) to al- 
low derotation to /I = 0. The 3cr error in position angle due to 
a Icr error of 1 rad m"^ at » 0.77 is 2.3 rad. This implies 
that the 99% confidence range is more than 4/3 times the 180° 
range that the position angle can assume. 

On the theoretical side we are faced with the following 
questions: 

- How can we explain the uniformity of the excess RM, in 
both sign and magnitude, across the face of the Perseus 
cluster? 

- How can we distinguish between the various models for 
the origin of the relativistic particles at the edges of the 
cluster (are they from re-energized relics or are they shock- 
accelerated?) 

- How long can AGN bubbles hibernate before they fail to be 
re-energized to radiate at observable radio frequencies? 

- What are the observable efifects if massive galaxies, or 
galaxy sub-groups, plough through such bubbles? The nu- 
merical simulations suggest that this would happen rather 
frequently. 

Bubbles carry important information on the final phases 
of radio galaxies and the study of AGN fossils may bring 
us somewhat closer to an answer on the long-standing ques- 
tion: do radio galaxies stop being active rather suddenly or 
do they just fade away? The study of bubbles and relics may 
quantify the role played by AGN activity and its feedback 
on cluster evolution. Cluster evolution and radio AGN activ- 
ity are clearly linked in the central parts of clus ters and may 
provide the answer to ' the cooling catastrophe' jpabian et alJ 
l2003H Clarke et al.'2004V Are enough bubbles blown to prevent 
such a catastrophe? It is be lieved that AGN activity is recurrent 
dSchoenmakers et alJ2000l) but it is at this moment unclear how 
many bubbles exist in a cluster, nor do we know how long they 
can hibernate. Is the large-scale polarized structure that arches 
around the steep spectrum tail of NGC 1265 indeed the remains 
of an earlier phase of activity? 

10. Conclusions and summary 

This paper has described low frequency observations of the 
Perseus cluster obtained with the WSRT. The total intensity 
data agree with earlier observations and do not present new in- 
sights into the cluster emission. However, using a novel method 



of wide band polarimetric imaging, called RM-synthesis, we 
have been able to extract information on extremely faint polar- 
ized signals buried deep in the noise of individual narrow band 
frequency images. 

Polarized emission has been detected in the Perseus clus- 
ter over a wide range of Faraday depths from about to 
90 rad m"'. Low emission ((f) < 15 rad m"^) is at- 
tributed to the local Galactic foreground. Emission at values 
of > 30 rad m"^ on the other hand, shows organized struc- 
tures on scales up to a degree and displays rapidly fluctuating 
polarization angles on scales of the order of l'-2'. A Galactic 
foreground interpretation for the high ^ emission appears ex- 
tremely implausible. The diffuse structures have a clear excess 
Faraday depth of about +40 rad m~^ relative to a dozen distant 
radio galaxies surrounding the Perseus cluster and are probably 
located at the near side periphery of the cluster. 

Most of the polarized emission, which is located at dis- 
tances of 0?5-l°5 from the cluster centre, appears 1-2 orders 
of magnitude too bright to be explainable as Thomson scattered 
emission o f the central radi o source off the thermal electrons in 
the cluster llSvunvaevll982l) . However, this remains a viable ex- 
planation for previously detected highly polarized 2 1 cm emis- 
sion from the inner 10'-20' as well as part of the 81-95 cm 
emission observed in the region surrounding NGC 1275. We 
believe that the bulk of the emission associated with the Perseus 
cluster is related to buoyant bubbles of relativistic plasma, 
probably relics from still active or now dormant AGN within 
the cluster. A lenticular shaped structure, referred to as the lens, 
and measuring 0. 5-1 Mpc is strikingly similar t o the struc- 
tures predicted by E nBlin & Gopal-Krish na ( 2001) and may be 
in the process of being re-energized. A bright, about 0.4 Mpc 
sized, structure located to the north of the steep spectrum tail of 
NGC 1265 may be evidence of an even earlier phase o f activity 
from this famous head-tail radio galaxy (iSiibring & de Bruvnl 
11998). 

We have not yet detected the total intensity corresponding 
to any of these polarized structures. This may be due to a com- 
bination of dynamic range in Stokes / and lack of short spacing 
sensitivity. However, there is no doubt that the structures are 
very highly polarized (>20%), as is common for many relic 
sources in clusters. 

At the western edge of the cluster we have detected lin- 
ear structures of several Mpc length that may be related to 
shocks caused by infall of gas into the Perseus cluster along 
the Perseus-Pisces filamentary structure of the cosmic web 
t.EnBlin et al. 1998: Burns 1998.) . 

The ability of low-frequency RM-synthesis to image very 
extended polarized structures, at intensity levels far below 
the total intensity image which may be limited by confusion 
or dynamic range problems, makes it a very powerful tech- 
nique. Regions of low surface brightness and low magneto- 
ionic plasma density are predicted to occur at the periphery of 
clusters and extending into the inter cluster medium and the in- 
terface with regions shocked in the process of large-scale struc- 
ture formation. The radio spectra associated with such emis- 
sion are often (predicted to be) steep requiring observations 
at low frequencies. The Low Frequency Array (LOFAR, see 
http://www.lofar.orgJ is the ideal instrument to embark on such 
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studies. It combines excellent surface brightness sensitivity and 
a wide frequency range with the capabiUty to perform high an- 
gular resolution studies in cases where beam depolarization is 
important. 
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